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We present the first computations of quasiequilibrium binary neutron stars with differ- 
ent mass components in general relativity, within the Isenberg- Wilson-Mathews approxima- 
tion. We consider both cases of synchronized rotation and irrotational motion. A poly- 
tropic equation of state is used with the adiabatic index 7 = 2. The computations have 
been performed for the following combinations of stars: (M/R)oo. star i vs. (A//i?) 00j sta r 2 = 
0.12 vs. (0.12, 0.13, 0.14), 0.14 vs. (0.14, 0.15, 0.16), 0.16 vs. (0.16, 0.17, 0.18), and 0.18 vs. 0.18, 
where (M/R)^ denotes the compactness parameter of infinitely separated stars of the same baryon 
number, ft is found that for identical mass binary systems there is no turning point of the binding 
energy (ADM mass) before the end point of the sequence (mass shedding point) in the irrotational 
case, while there is one before the end point of the sequence (contact point) in the synchronized 
case. On the other hand, in the different mass case, the sequence ends by the tidal disruption of 
the less massive star (mass shedding point). It is then more difficult to find a turning point in 
the ADM mass. Furthermore, we find that the deformation of each star depends mainly on the 
orbital separation and the mass ratio and very weakly on its compactness. On the other side, the 
decrease of the central energy density depends on the compactness of the star and not on that of 
the companion. 

PACS numbers: 04.25.Dm, 04.40.Dg, 97.60.Jd, 97.80.-d 



I. INTRODUCTION 



Coalescing binary neutron stars are expected to be one of the most promising sources of gravitational waves 
that could be detected by the ground based, kilometer size laser interferometers such as the Laser Interferometric 
Gravitational wave Observatory (LIGO), VIRGO, GEO600, and TAMA300 1 , and also are considered as one of 
candidates of gamma-ray burst source . 

Due to the emission of gravitational radiation, binary neutron stars decrease their orbital separations and finally 
merge. When discussing the evolution of the system, it is convenient to separate it into three phases. The first 
one is the inspiraling phase in which the orbital separation is much larger than the neutron star radius, and the 
post-Newtonian (PN) expansion constitutes an excellent approximation. Recently, two groups succeeded in deriving 
the 3PN equation of motion of point-mass binary systems j}[ f| , and the equivalence of the results between these 
groups is shown in The second stage is the intermediate phase in which the orbital separation becomes only 

a few times larger than the radius of a neutron star, so that hydrodynamics as well as general relativity play an 
important role. In this phase, since the shrinking time of the orbital radius due to the emission of gravitational waves 
is still larger than the orbital period, it is possible to approximate the state as quasiequilibrium. The final stage is 
the merging phase in which the two stars coalesce dynamically. As in the intermediate phase, since hydrodynamics 
as well as general relativity play an important role, fully relativistic hydrodynamical treatments are required in this 
phase which therefore pertains to the field of numerical relativity The first successful computations of the evolution 
of binary neutron stars from their innermost stable circular orbits (ISCO) to black hole or massive neutron star 
formation have been performed by Shibata and Uryu [R | |. Other efforts are presented in [|[ |Io| . 

The present article belongs to the series of works fniril (hereafter Paper I), Jl|| (hereafter Paper II), |J], devoted 
to the intermediate phase. This stage is interesting because we may get informations about equation of state of neutron 
stars through the ISCO |15| ]. Furthermore, it is important from a numerical point of view since it provides initial data 
for the dynamical simulation in the merging phase |t], |^] . Then numerous theoretical efforts are devoted in this phase 



1 For TAMA300, data taking has started in summer 1999, and the first results of the data analysis have been published [ 
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including (semi-) analytic Newtonian Jh| [l7| and post-Newtonian [fill |l9l , |20[ ap proaches and numerical Newtonian 
H, H, |2§ ||, post-Newtonian and general relativistic fnj/pk g§ |27J §|, ||, |l], || ones. Among these 
studies, we concentrate on numerical computations in the general relativistic framework. Kochaneck fl33| and Bildsten 
and Cutler |]34f have shown that the gravitational-radiation driven evolution is too rapid for the viscous forces to 
synchronize the spin of each neutron star with the orbit as they do for ordinary stellar binaries. Rather, the viscosity 
is negligible and the fluid velocity circulation (with respect to some inertial frame) is conserved in these systems. 
Provided that the initial spins are not in the millisecond regime, this means that close binary configurations are well 
approximated by zero vorticity (i.e. irrotational) states. Formulations for irrotational binaries in general relativity 
have been developed by several authors |35], [57| [38| (see Appendix A of Paper 1 for a comparison between them) . 
In addition to the quasiequilibrium and irrotational assumptions, we use the approximation of a conformally flat 
spatial 3-metric, introduced by Isenberg |39| and Wilson and Mathews Q] (hereafter IWM approximation; see Sec. 
IV. C of Q] and Sec. III. A of Paper I for a discussion). 

Until now, several groups have produced quasiequilibrium configurations of binary neutron stars, as listed above. 
Among them, there are results for synchronized |^6|, |2?J and irrotational jll], [l^, |||, |9[ rotation states in the general 
relativistic framework, within the IWM approximation, and those for synchronized one within some axisymmetric 
approximation |3l[ |32| . However they all deal with identical star binaries and calculations for different mass binary 
system in general relativity have not been performed yet (in Newtonian theory they have been performed in the 
synchronized case |l4|, |2j|, and recently in the irrotational one [p^[). 

We will present here the first numerical results of relativistic binary systems composed of different mass stars. In 
addition, we give results for systems of identical stars which extend those already presented in Paper I. As a first 
step in the study of different mass systems, we restrict our computations to a polytropic equation of state, with the 
adiabatic index 7 = 2. 

The plan of the article is as follows. A brief overview of our method is given in 

Sec.H 

In Sec. IIR new tests of the 



numerical code are shown, which were not given in Papers I and II. We present the numerical results in Sec. IV, and 
discuss them in Sec. [v|. Finally Sec. VI is devoted to the summary. 



Throughout this article, we adopt geometrical units: G = c = 1, where G and c denote the gravitational constant 
and speed of light, respectively. 

II. METHOD 

The reader is referred to Sees. II, III, and IV of Paper I for the complete set of equations governing perfect fluid 
binary stars under the assumption of quasiequilibrium and conformally flat spatial metric, as well as the numerical 
method which we use. We simply give here some outline of the method. 

Let us summarize the assumptions expressed in Sec. [|. The first one is quasiequilibrium, resulting from the fact 
that the time scale of orbital shrinking is larger than that of the orbital revolution until the ISCO. This assumption 
is taken into account by demanding that the spacetime is endowed with a helical Killing vector 

where f2 is the orbital angular velocity and the vectors d/dt and 8/ dip are associated with the time coordinate t 
and the azimuthal coordinate (p of an asymptotical inertial observer. The second assumption regards the matter 
stress-energy tensor, which we assume to have the perfect fluid form: 

Tpv = (e + p)u^u u + pg^, (2) 

where e denotes the fluid proper energy density, p the fluid pressure, the fluid 4-velocity, and g^ the spacetime 
metric. The third assumption concerns the rotation state of the binary system. Although the realistic rotation state 
will be an irrotational one (cf. Sec. |), we consider both cases of synchronized and irrotational states in order to 
exhibit their differences. The fourth assumption is on the equation of state. For simplicity, we use a polytropic one, 

p = nn 1 , (3) 

where n is the fluid baryon number density and k and 7 are some constants. Here we assume that the two neutron 
stars have the same equation of state, i.e., the constants K and 7 are identical for both stars. The fifth and last 
assumption is the IWM one, namely a conformally flat spatial metric. Then the full spacetime metric takes the form: 

ds 2 = -(A 2 - B.B^dt 2 - 2B i dtdx i + ,l-/, ,(/r'-/.r'\ (4) 

N being the lapse function, B l the shift vector, A the conformal factor, and f^ the flat spatial metric. Under these 
assumptions, the equations for fluid motion and gravitational field are presented in Paper I. 
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FIG. 1: Relative error in the virial theorem for numerical solutions of relativistic quasiequilibrium binary system, as a function of 
the orbital separation. Solid, dashed, and dotted lines denote the relative errors defined by Eq. (0), (0) and (|l|), respectively. 
Left panel is for synchronized binaries, and the right one for irrotational ones. Both of them are drawn for identical mass binary 
systems with the compactness M/R = 0.14 vs. 0.14. 



III. CODE TESTS 

In this series of studies jy], [l3| Q, we use a multidomain spectral method with surface-fitted coordinates 
p2| , |43| , ft"4| . This method is implemented by a numerical code constructed upon the C++ library Lorene p5| . 
Numerous tests of the code have been already presented in Papers I and II. In particular, the comparison with results 
from other groups j2(| ^9) has been given in Sec. V.D of Paper I for the identical mass case. We will give here new 
tests based on the relative error in a general relativistic generalization of the virial theorem. 

In Newtonian gravity, the virial theorem has proved to be useful to check the global error in numerical solutions for 
stationary fluid systems. A general relativistic version of the virial theorem has been obtained by Gourgoulhon and 
Bonazzola jl6| for stationary spacetimes. It has been recently extended to binary star spacetimes within the IWM 
approximation by Friedman, Uryu and Shibata [^lj 2 . The virial relation is equivalent to 

-MADM - Alomar = 0, (5) 

where Madm is the ADM mass: 

A^adm = j V l A 1 / 2 dS l , (6) 
and MKomar is a Komar-type mass, defined by 

AfKomar = ^ j ^NdS*. (7) 

The virial theorem obtained by Friedman, Uryu and Shibata |4l| (see also Eq. (5.7) of Ref. |47]]) writes 

VE{FUS) = j \2NA 3 S + -^-NA 3 K 3 i K i j + l~NA<^ t 0Tp - V^VV) 
J I on An 

= 0. (8) 



2 The spacetime generated by a binary system is not stationary, due to gravitational radiation. However, within the IWM approximation, 
the gravitational radiation is neglected in the global spacetime dynamics, so that one is able to recover the virial theorem. 
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By a straightforward manipulation, this integral can be recast in the form of the virial theorem as obtained by 
Gourgoulhon and Bonazzola |46[: 



VE(GB) = 



2A 3 S + it-J^KiK) + -^A(V l (3V l l3 - V^VV - 2V./3VV) 

07T J 47T 



0. 



(9) 



Let us stress that the above identity has been derived by Gourgoulhon and Bonazzola j46j only for stationary space- 
times and that its validity for IWM spacetimes with helical Killing vector has been obtained by Friedman, Uryu and 
Shibata 

As error indicators of our numerical solutions, we have evaluated the quantities 



Mydm - M K O 



DM 



(10) 
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The results are presented in Fig. [j] as a function of the orbital separation. One can see from this fi; 
typical virial error is ~ 1 x 10 -5 through the sequence. 
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FIG. 2: Relative change in central energy density as a function of the orbital separation between the centers of mass of each 
star. Left (resp. right) panel is for synchronized (resp. irrotational) binaries with an adiabatic index 7 = 2. Solid line denotes 
the identical mass case M/R = 0.14 vs. 0.14. The results of the different mass case M/R = 0.14 vs. 0.16 are shown as dotted 
(less massive star) and dashed (massive star) lines. 



IV. NUMERICAL RESULTS 



In this section, we present the numerical results for quasicquilibrium sequences of constant baryon number binaries 
(evolutionary sequences), for both synchronized and irrotational cases and for the adiabatic index 7 = 2. Although 
the interesting case is the irrotational one, as mentioned in Sec. [j], we have computed configurations for both rotation 
states in order to compare their properties. The number of computational domains is 4 or 5 for each star and the 
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space around it. The domains have an onion-like structure, being centered on each star: the innermost domain covers 
the star, two or three domains (shells) are placed outside it, and the external domain extends to infinity thanks to 
some compactification (cf. Sec. IV. A of Paper I for all details). The total number of domains is therefore 2 x 4 = 8 or 
2 x 5 = 10. We use 8 domains for close orbital separation and 10 domains for large orbital separation. The numbers 
of collocation points in each domains are N r x Ng x N v — 33 x 25 x 24 or 25 x 17 x 16, where N r , Ng, and N v denote 
the number of collocation points of the radial, the polar, and the azimuthal directions, respectively. 

We parametrize the different mass systems by the values of (M/R) oa> stal - i vs. (M/R)^. sta r 2, where (M/R)^ 
denotes the compactness of the isolated spherical star with the same baryon number. 

First of all, the relative changes in central energy density for the cases of M/R = 0.14 vs. 0.14 and 0.14 vs. 0.16 are 
depicted in Fig. ||as a function of the orbital separation. (Hereafter we abbreviate (M/R) 00< star i vs. (M/R) 00t sta r 2 
as M/R for simplicity.) One can see from these figures that the central energy density decreases monotonically in 
both synchronized and irrotational cases 3 . The decrease of the central energy density for the synchronized case is 
about 10 % or more at the end of the sequence, while for the irrotational case, it is only about 1 %. 
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FIG. 3: ADM mass of the binary system as a function of the orbital angular velocity. Left (resp. right) panel is for synchronized 
(resp. irrotational) binaries. In each panel, the three sub-panels correspond to the compactness of M/R = 0.12 vs. 0.12, 0.12 
vs. 0.13, and 0.12 vs. 0.14 from the left to the right. 




3 In Paper I (Table III), as well as in Jll| (Fig. 1), there was a slight increase in the central energy density for irrotational binaries at the 
large and medium separations. However, from Paper II we have improved the method to determine the surface of the star, which leads 
to a better accuracy in the present work. 



FIG. 4: Same as Fig. |, but for the compactness of M/R = 0.14 vs. 0.14, 0.14 vs. 0.15, and 0.14 vs. 0.16 from the left to the 
right in each panel. 




FIG. 5: Same as Fig. |, but for the compactness of M/R = 0.16 vs. 0.16, 0.16 vs. 0.17, and 0.16 vs. 0.18 from the left to the 
right in each panel. 

Next, the ADM mass of the binary system [Eq. (||)] is shown in Figs. ||- [s]as a function of the orbital angular velocity, 
for different mass ratios and different states of rotation. One can clearly see the turning point in the synchronized 
cases for identical mass binaries. For more compact stars the turning points appear relatively earlier than those of 
less compact ones. On the other hand, we do not find any turning point in irrotational binaries. Moreover, it becomes 
more difficult for the turning point to appear even in the synchronized case when the mass difference in the binary 
system increases. 

In Fig. ^, we show the total angular momentum of the binary system as a function of the orbital angular velocity for 
the compactness M/R = 0.18 vs. 0.18. For synchronized binaries, the turning point of the total angular momentum 
coincides with that of the ADM mass within the numerical accuracy. 
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FIG. 6: Total angular momentum of the binary system as a function of the orbital angular velocity for the compactness 
M/R = 0.18 vs. 0.18. Left panel is for the synchronized case and right for the irrotational one. 
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FIG. 7: Equatorial to polar ratio of the radial derivative of the enthalpy as a function of the orbital separation. Left (resp. 
right) panel is for synchronized (resp. irrotational) binaries. Solid, dotted, dashed, and long-dashed lines denote the cases of 
the compactness M/R = 0.12 vs. 0.12, 0.14 vs. 0.14, 0.16 vs. 0.16, and 0.18 vs. 0.18, respectively. 
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FIG. 8: Same as Fig. [j], but for different mass binaries. Left (resp. right) panel is for synchronized (resp. irrotational) binaries. 
Solid line denotes the case of identical mass binary with the compactness M/R = 0.14. Thin and thick dotted lines are for the 
compactness M/R — 0.14 vs. 0.15 of the more massive and less massive stars, respectively. Thin and thick dashed lines are for 
the compactness M/R — 0.14 vs. 0.16 of the more massive and less massive stars, respectively. 



The equatorial to polar ratios of the radial derivative of the enthalpy are shown as a function of the orbital separation 
normalized by the radius of a star to the companion star in Figs, [z] and |[ As explained in Papers I and II, this 
quantity 



X ■= 



(dH/dr) 



cq,comp 



(dH/dr) pole 



(13) 



is a good dimensionless indicator to detect the appearance of cusp at the stellar surface, which corresponds to the 
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mass shedding limit. It is worth to plot the variation of \ with respect to the orbital separation normalized by the 
sum of the radius of each star toward its companion: dj (ai + a'^J, where d denotes the orbital separation between the 
two coordinate center (maximum of density), a\ the radius of the less massive star in the direction of its companion, 
and a[ that of the more massive star. Note here that the separation d and the radii a\ and a[ are coordinate lengths. 
If the quantity x becomes zero, we can conclude that the sequence ends by mass shedding. On the other hand, if the 
quantity d/{a\ + a'i) becomes unity, the two stars will contact each other. It is found from Fig. [7] that for synchronized 
identical mass binaries the sequences seem to terminate by the contact between the two stars (the quantity x will 
become zero at the same time). The important feature of the sequences is that they take almost the same track 
regardless of the compactness. On the contrary, one can see from Fig. [?] that for irrotational binaries the sequences 
seem to terminate by a cusp point (mass shedding). Such a behavior is enhanced for more compact cases. In Fig. ||, we 
compare the quantity x of identical mass binary of the compactness M/R = 0.14 with that of different mass binaries 
of the compactness M/R = 0.14 vs. 0.15 and 0.14 vs. 0.16. It is found that the sequences seem to terminate by a 
cusp point of the less massive star for both synchronized and irrotational binaries when we extrapolate the lines, while 
the more massive star keeps its shape near spherical one (see the next paragraph). When we increase the difference 
in mass between two stars, this tendency is enhanced. 

It is worth to mention that as soon as the two masses differ, even slightly, the sequence of a synchronized binary 
system terminates by a cusp point, because the line for the less massive star is always located below that for the 
identical mass stars which will reach x — at d/ (a\ + a[) = 1. This implies that the line for the less massive star 
seems to always reach x = at d/(a\ + a[) > 1. Note here that the sequences for identical mass binaries and less 
massive stars in different mass binaries should reach the point x — 0- However, in our numerical method, we cannot 
treat a cuspy figure, because we adapt the innermost numerical domain to the surface of the star and the multidomain 
spectral method that we use assume that the domain boundaries are smooth surfaces. For very close configurations, 
it becomes therefore difficult to make the iterative procedure converge to a sufficient level, i.e. to have the relative 
difference SH in the enthalpy field between two steps lower than 10 -5 . This explains why we stop the sequences at 
X — 0.4. We use SH = 10~ 10 for very large separations, SH = 10 -8 — 10~ 7 for medium ones, and SH = 10~ 6 — 10~ a 
for very close ones. 



Baryon density (y=0) Baryon density (y=0) Baryon density (y=0) 




FIG. 9: Isocontours of the baryon density for synchronized binary systems. Left panel is for the compactness M/R = 0.14 vs. 
0.14, the center panel M/R = 0.14 vs. 0.15, and the right panel M/R — 0.14 vs. 0.16. The orbital axis is located on x = in 
each panel. 

In Figs. |^ and [To], isocontours of the baryon density of various binary systems are shown. There are three panels 
in each figure. The left panel shows the result for identical mass binary with the compactness M/R — 0.14 vs. 0.14. 
The center and right ones show the results for different mass binaries with the compactness M/R = 0.14 vs. 0.15, and 
0.14 vs. 0.16, respectively. In these panels, the left hand side star is the less massive star. It is clearly seen that the 
less massive star is tidally deformed and elongated while the more massive star relatively does not deviate from the 
spherical shape so much (see the right panel) . The orbital separations of these isocontours are listed at the last line of 
each case in Tables below which corresponds to the closest figure we can calculate, except for the case ofM/i? = 0.14 
vs. 0.15 which we do not show in the present paper as a table. 

Finally, our numerical results for constant baryon number sequences are presented in Tables || - IV . We represent 
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FIG. 10: Same as Fig. bl but for irrotational binary systems. 



the orbital separation by two quantities: 



do = 
d = 



dg 

-^poly 
d 



R 



poly 



(14) 
(15) 



da is the coordinate separation between the centers of mass of each star (same definition as Eq. (107) of |^9| or Eq. 
(128) of Paper I), d that between the two stellar centers defined as the maxima of the density field, and i? po i y the 
length constructed from the polytropic constants k and 7 [Eq. (0)] 



-^poly ■ — K ^ ^ ^ . 

The dimensionless orbital angular velocity is defined by 

Q:= $7^/2(7-1), 

the dimensionless ADM mass by 

M := AfADM^ 172 ^- 1 ', 



the dimensionless baryon mass by 



Mr := MbK- 1 /*" 1 ' 



and the dimensionless total angular momentum by 

J:= Jk- 1 ^- 1 ). 



(16) 
(17) 
(18) 
(19) 
(20) 



The virial errors VE(M), VE(GB) and VE(FUS) are defined by Eqs. @, |[J) and (|TT|) respectively. a x denotes 
the coordinate radius toward the companion star, 02 the radius perpendicular to it in the orbital plane and CZ3 the 
radius perpendicular to the orbital plane. Se c is the relative change in central energy density: Se c :— (e c — ecooVecoo, 
and the cusp indicator x is defined by Eq. ([l3]). In the Tables, the superscript " ' " denotes the values for the more 
massive star, and the symbol t denotes the turning point in the ADM mass along the sequence. 



V. DISCUSSION 



In this section, we discuss about the behaviors of the axial ratios, the relative change in central energy density, and 
the turning point of the synchronized sequences. 
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FIG. 11: Axial ratios 02/01 or 03/01 as a function of the cusp indicator \. Left (resp. right) panel is for synchronized (resp. 
irrotational) binaries. Solid and dotted lines denote the axial ratios 02/01 and 03/01 for the compactness M/R — 0.12 vs. 0.12, 
respectively. Dashed and long-dashed lines are the axial ratios 02/01 and 03/ai for the compactness M/R — 0.16 vs. 0.16, 
respectively. 



In Figs. |Tl| and [12|, the axial ratios are shown as a function of the cusp indicator \- I n Fig- 11 , we show the axial 
ratios 02/01 and 03/01 for identical mass binaries with the compactness M/R = 0.12 vs. 0.12 and 0.16 vs. 0.16. One 
can see from this figure that each axial ratio has almost the same track regardless of the difference of the compactness, 
in particular for synchronized binary systems. Moreover, such a coincidence of the sequence appears in the cases of 
different mass binary systems (see Fig. |l2]). In those cases, the sequences almost coincide with each other regardless 
of the compactness of the companion star. Note that the sequences of the axial ratios for more massive stars end 
before x — 0, while those for less massive stars should reach the point x = 0. However we cannot treat a cuspy figure 
with our numerical method, because we adapt the innermost numerical domain to the surface of the star and the 
domain boundaries are assumed to be smooth (see discussion in Sec. IV). Therefore we had to stop the sequences at 
around \ — 0-4 for the less massive star. 

The relative change in central energy density is presented as a function of \ m Fig- It is clearly seen that the 
sequences for the synchronized case with the same compactness coincide with each other, while those with the different 
compactness split into curves with respect to each compactness. For example, the relative change for the identical 
binary system with the compactness M/R = 0.14 coincides with that for the more massive star of the different mass 
binary with M/R = 0.12 vs. 0.14 and that for the less massive star with M/R = 0.14 vs. 0.16. However, the 
sequences for stars with M/R = 0.12, 0.14, and 0.16 are different from each other. 

Finally, let us comment on the behavior of the turning point of the ADM mass (binding energy) along a sequence. 
We never found any turning points for irrotational binary systems with the polytropic index 7 = 2. On the other hand, 
it appeared clearly for the synchronized case. In Fig. 14, the cusp indicator \ at the turning point of synchronized 
identical mass binary systems is plotted as a function of the compactness M/R. Since the accuracy of the determination 
of the ADM mass is order of 10 -5 in our calculations, the relative errors on the orbital angular velocity and the quantity 
X at the turning point become a few percents. Even when one is taking this fact into account, it is possible to conclude 
that the line in Fig. |TJ increases proportionally to the compactness of the stars. This behavior agrees with the value 
of Xtuming that we have found in Newtonian calculations (M/R — 0), namely Xturning = 0.2862 (Table I of Paper II). 



VI. SUMMARY 



In the present paper, we have computed constant baryon number sequences of binary neutron stars in quasiequilib- 
rium, in both cases of synchronized and irrotational motion and in both cases of identical and different mass systems. 
We have performed a general relativistic treatment, within the IWM approximation (conformally flat spatial metric). 



11 



Axial ratios 

Synchronized case, y=2 



0.9 



0.7 



0.6 



IiTolalional case, y=2 





-j M/R=0 
> M/R=0 
□ M/R=0 
■ M/R=0 
j M/R=0 
» M/R=0 
iM/R=0 
i M/R=0 
-7 M/R=0 
f M/R=0 



14vs.0.14: 
14vs.0.14: 
14vs.0.15: 
14vs.0.15: 
14vs.0.15: 
14vs.0.15: 
14vs.0.16: 
14vs.0.16: 
14vs.0.16: 
14vs.0.16: 



(0.14) 
(0.14) 
(0.15) 
(0.15) 
(0.14) 
(0.14) 
(0.16) 
(0.16) 



0.2 



0.4 0.6 

X 



0.2 



0.4 0.6 0.8 1 

X 



FIG. 12: Same as Fig. [llj but for different mass binaries. Solid line with open and filled circles denote the axial ratios ai/<X\ and 
(23/0,1 of the compactness M/R = 0.14 vs. 0.14, respectively. Dotted line with open and filled squares denote the axial ratios 
02/01 and 03/01 of the compactness M/R = 0.14 vs. 0.15 for less massive star, respectively. Dashed line with open and filled 
diamonds denote the axial ratios 02/01 and 03/01 of the compactness M/R = 0.14 vs. 0.15 for more massive star, respectively. 
Long-dashed line with open and filled triangle- up denote the axial ratios 02/01 and 03/01 of the compactness M/R = 0.14 vs. 
0.16 for less massive star, respectively. Dot-dashed line with open and filled triangle-down denote the axial ratios 02/01 and 
03/01 of the compactness M/R = 0.14 vs. 0.16 for more massive star, respectively. 



We have used a polytropic equation of state with an adiabatic index 7 = 2. 

The summary of our results is that (1) Among the quasiequilibrium sequences we have calculated in the present 
article, only the synchronized identical star binaries terminate by the contact between two stars, while all the other 
types of sequences (synchronized different star binaries, irrotational identical star and different star binaries) end at 
the mass shedding points. Note here that we cannot conclude about the end points of non- quasiequilibrium sequences 
(i.e. dynamical sequences), such that binary systems without viscosity but with intrinsic non-aligned spins, those 
with slight viscosity and deviating from both the irrotational state and the synchronized one, those with infalling 
radial velocity, and so on. (2) For identical mass binary systems there is no turning point of the ADM mass in the 
irrotational case, while there is clearly one before the end point of the sequence in the synchronized case. (3) It is 
more difficult to see the turning points of the ADM mass for different mass binary systems than for identical ones. 
(4) The deformation of the star is determined by the orbital separation and the mass ratio and is not affected much 
by its compactness. (5) The decrease of the central energy density depends on the compactness of the star and not 
on that of its companion. 
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TABLE I: Orbital angular velocity, ADM mass, total angular momentum, virial errors, axial ratios, relative change in central 
energy density, and the cusp indicator \ along constant baryon number sequences of identical mass binaries in synchronized 
motion, for various compactness. The symbol t denotes the turning point of the sequences. 
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3100 


9.567{ 


-2) 




4.965 


-5) 


1.943 


-4) 


1.632 


-4) 


0.8537 





7930 


-1 


747( 


-1) 





4120 
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TABLE II: Same as Table | but for different mass binary systems. The superscript " ' " denotes the values for the more massive 
star. 



do 


d 


Q 




M 




Different mass stars, Synchronized case 
J \VE~{M)\ \VE{GB)\ \VE(FUS)\ 


Q>2 1 Q-i 

O,o / (2i 


0,3 J d\ 
i i f 




Sc c 






A. 
v' 

A 














M/R = 


0.12 vs 


0.14, 


M B 


= 0.1299 vs. 


1461 
















4.873 


2.863 


4.404 


-2) 


0.2570 


8 


347(-2) 


1.201( 


-5) 


2.022 


-5) 


1.895( 


-5) 


0.9827 


0.9754 


_1 


197 


_2) 





9345 




























0.9891 


9833 


_1 


066 


_2) 


o 


9548 


3.655 


2.117 


6.621( 


-21 


0.2567 


7 


675 (-2) 


9.690( 


-6) 


1.633 


-5) 


1.501( 


-5) 


9593 


0.9391 


_2 


824 


_2) 


o 


8376 




























0.9750 


9592 


_2 


467 


_2) 


o 


8891 


3.046 


1.715 


8.547 


-21 


2565 


7 


382(-2) 


7.292( 


-6) 


1.220 


-5) 


1.106( 


-5) 


0.9234 


8882 


_5 


023 


_2) 


o 


7001 




























0.9547 


9265 


_4 


264 


_2) 


o 


8001 


2.864 


1.583 


9.317 


-2) 


0.2564 


7 


312(-2) 


6.226( 


-6) 


1.043 


-5) 


9.430( 


-6) 


0.9018 


0.8600 


-6 


189 


_2) 


o 


6227 




























0.9434 


9095 


_5 


169 


_2) 


o 


7532 


2.742 


1.488 


9.899 


-2) 


0.2564 


7 


275(-2) 


1.408( 


-5) 


1.812 


-5) 


1.523( 


-5) 


0.8807 


0.8338 


_7 


209 


_2) 


o 


5501 




























9330 


0.8946 


_5 


930 


_2) 


o 


7121 


2.681 


1.437 


1.022 


-1) 


0.2563 


7 


261(-2) 


1.605( 


-5) 


2.080 


-5) 


1.631( 


-5) 


0.8665 


0.8170 


_7 


840 


_2) 


() 


5026 




























0.9265 


0.8856 


-6 


392 


_2) 


o 


6870 


2.620 


1.383 


1.055 


-1) 


2563 


7 


249(-2) 


1.185( 


-5) 


2.514 


-5) 


2.050( 


-5) 


0.8486 


0.7965 


-8 


565 


_2) 


o 


4439 




























0.9190 


0.8753 


-6 


900 


_2) 


o 


6582 


2.577 


1.342 


1.080 


-I) 


2563 


7 


243(-2) 


1.161( 


-6) 


4.051 


-5) 


3.446( 


-5) 


0.8327 


0.7789 


_g 


144 


_2) 


o 


3923 




























0.9128 


0.8671 


_7 


294 


_2) 


o 


6352 














M/R = 


0.14 vs 


0.16, 


Mb 


= 0.1461 vs. 


1600 
















4.873 


3.102 


4.586 


-2) 


0.2823 


9 


688(-2) 


2.076( 


-5) 


3.649 


-5) 


3.473( 


-5) 


0.9865 


0.9810 


-1 


164 


_2) 


o 


9484 




























0.9912 


0.9869 


-1 


119 


_2) 


o 


9635 


3.655 


2.307 


6.879 


-2) 


0.2818 


8 


905 (-2) 


1.736( 


-5) 


3.055 


-5) 


2.855( 


-5) 


0.9694 


0.9534 


_2 


691 


_2) 


() 


8733 




























9805 


0.9679 


_2 


550 


_2) 


o 


9110 


3.046 


1.888 


8.858 


-2) 


0.2816 


8 


552(-2) 


1.370( 


-5) 


2.420 


-5) 


2.232( 


-5) 


0.9449 


0.9161 


_4 


664 


_2) 


o 


7717 




























0.9657 


0.9429 


_4 


331 


_2) 


o 


8417 


2.742 


1.662 


1.024 


-1) 


0.2814 


8 


408(-2) 


1.132( 


-5) 


1.941 


-5) 


1.771( 


-5) 


0.9186 


0.8794 


-6 


518 


_2) 


o 


6711 




























0.9507 


0.9193 


_5 


931 


_2) 


o 


7762 


2.620 


1.564 


1.090 


-1) 


0.2814 


8 


364(-2) 


2.342( 


-5) 


3.453 


-5) 


2.770( 


-5) 


0.9016 


0.8572 


_7 


584 


_2) 


o 


6096 




























0.9416 


9056 


-6 


817 


_2) 


o 


7384 


2.499 


1.460 


1.163 


-I) 


0.2814 


8 


329(-2) 


2.577( 


-5) 


3.575 


-5) 


2.861( 


-5) 


0.8774 


0.8271 


-8 


945 


_2) 


() 


5251 




























0.9295 


0.8882 


_7 


911 


_2) 


o 


6897 


2.438 


1.403 


1.203 


-I) 


0.2813 


8 


318(-2) 


1.268( 


-5) 


4.801 


-5) 


3.990( 


-5) 


8603 


0.8071 


_g 


800 


_2) 


o 


4677 




























0.9217 


0.8775 


-8 


567 


_2) 


o 


6595 


2.377 


1.341 


1.245 




0.2813 


8 


310(-2) 


1.713( 


-5) 


1.014 


-4) 


8.814( 


-5) 


0.8376 


0.7814 


_1 


079 


_1) 


o 


3927 




























0.9124 


8650 


_g 


308 


_2) 


o 


6238 














M/R = 


0.16 vs 


0.18, 


Mb 


= 0.1600 vs. 


1710 
















4.264 


2.957 


5.704 


-2) 


0.3023 


1 


037(-l) 


3.359( 


-5) 


6.089 


-5) 


5.824( 


-5) 


0.9852 


0.9782 


-1 


745 


_2) 


o 


9393 




























0.9901 


0.9848 


_1 


904 


_2) 


o 


9566 


3.655 


2.526 


7.074 


-2) 


0.3020 


9 


936(-2) 


3.049( 


-5) 


5.543 


-5) 


5.256( 


-5) 


0.9772 


0.9648 


_2 


712 


_2) 


o 


9021 




























9850 


0.9754 


_2 


934 


_2) 


o 


9301 


3.046 


2.080 


9.094 


-2) 


0.3017 


9 


529(-2) 


2.522( 


-5) 


4.605 


-5) 


4.311( 


-5) 


0.9602 


0.9373 


_4 


608 


_2) 


() 


8260 




























0.9743 


0.9564 


_4 


896 


_2) 


o 


8764 


2.620 


1.747 


1.116 


-1) 

*-} 


0.3015 


9 


293(-2l 


1.940( 


-5) 


3.529 


-5) 


3.267( 


-51 


0.9328 


0.8963 


_7 


264 


_2) 


o 


7127 




























0.9577 


9289 


_7 


522 


_2) 


o 


7990 


2.438 


1.592 


1.230 


-I) 


0.3014 


9 


221(-2) 


3.520( 


-5) 


5.539 


-5) 


4.539( 


-5) 


0.9100 


8653 


_g 


132 


_2) 


o 


6266 




























0.9448 


0.9089 


-9 


288 


-2) 





7431 


2.316 


1.479 


1.318 


-1) 


0.3013 


9 


188(-2) 


3.541( 


-5) 


5.009 


-5) 


4.242( 


-5) 


0.8862 


0.8350 


-1 


085 


-1) 





5411 




























0.9324 


0.8907 


-1 


081 


-1) 





6915 


2.237 


1.398 


1.381 


-1) 


0.3013 


9 


177(-2) 


1.418( 


-5) 


7.574 


-5) 


6.303( 


-5) 


0.8633 


0.8075 


-1 


226 


-1) 





4625 




























0.9214 


0.8753 


-1 


204 


-1) 





6480 


2.194 


1.350 


1.417 


-1) 


0.3013 


9 


174(-2) 


1.249( 


-5) 


1.266 


-4) 


1.084( 


-4) 


0.8465 


0.7883 


-1 


315 


-1) 





4063 




























0.9141 


0.8655 


-1 


278 


-1) 





6200 
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TABLE III: Same as Table jl but for irrotational binaries. 



Identical mass stars, Irrotational case 







d 




Q, 






M 




J 


\VE(M)\ 


VE(GB) 


\VE(FUS)\ 


aa/oi 






5e c 






X 




















M/R = 


0.12 vs 


0.12, 


Mb 


= 


1299 vs. 


1299 




















6.090 


3 


472 


3 


120( 


-2) 





2434 


8 


076(-2) 


9.264( 


-6) 


1.530 


-5) 


1.439( 


-5) 





9916 





9945 


-8 


061 


-5) 





9837 


4.872 


2 


780 


4 


304( 


-2) 





2432 


7 


365(-2) 


8.394( 


-6) 


1.371 


-5) 


1.266( 


-5) 





9846 





9882 


-2 


038 


-4) 





9654 


3.653 


2 


072 


6 


488 ( 


-2) 





2428 


6 


589(-2) 


5.992( 


-6) 


9.535 


-6) 


8.445( 


-6) 





9631 





9681 


-8 


804 


-4) 





9055 


3.044 


1 


694 


8 


397 


-2) 





2425 


6 


180(-2) 


3.377( 


-6) 


4.933 


-6) 


3.962( 


-6) 





9299 





9369 


-2 


613 


-3) 





8082 


2.861 


1 


570 


9 


164 


-2) 





2423 


6 


058(-2) 


6.164( 


-6) 


8.369 


-6) 


6.394( 


-6) 





9100 





9183 


-3 


929 


-3) 





7459 


2.738 


1 


478 


9 


746( 


-2) 





2423 


5 


978(-2) 


4.595( 


-6) 


6.938 


-6) 


5.079( 


-6) 





8896 





8988 


-5 


328 


-3) 





6692 


2.616 


1 


375 


1 


040( 


-1) 





2422 


5 


900(-2) 


1.337( 


-5) 


1.239 


-5) 


7.816( 


-6) 





8582 





8705 


-7 


548 


-3) 





5655 


2.554 


1 


304 


1 


076 ( 


-1) 





2422 


5 


864(-2) 


5.612( 


-7) 


1.134 


-6) 


6.392( 


-6) 





8283 





8428 


-9 


566 


-3) 





4476 


2.523 


1 


259 


1 


095( 


-1) 





2421 


5 


846(-2) 


2.005( 


-5) 


5.038 


-6) 


1.880( 


-5) 





8064 





8223 


-1 


111 


-2) 





3720 




















M/R = 


0.14 vs 


0.14, 


Mb 


= 


1461 vs. 


1461 




















4.872 


3 


001 


4 


510( 


-2) 





2707 


8 


755(-2) 


1.494( 


-5) 


2.565 


-5) 


2.409( 


-5) 





9877 





9914 


-1 


870 


-4) 





9736 


3.653 


2 


249 


6 


780( 


-2) 





2702 


7 


858(-2) 


1.136( 


-5) 


1.924 


-5) 


1.748( 


-5) 





9717 





9766 


-7 


114 


-4) 





9284 


3.044 


1 


856 


8 


753( 


-2) 





2699 


7 


384(-2) 


7.459( 


-6) 


1.222 


-5) 


1.054( 


-5) 





9483 





9545 


-1 


940 


-3) 





8594 


2.739 


1 


645 


1 


014( 


-1) 





2697 


7 


146(-2) 


5.141( 


-6) 


6.971 


-6) 


5.511( 


-6) 





9231 





9309 


-3 


673 


-3) 





7812 


2.617 


1 


553 


1 


080( 


-1) 





2696 


7 


053(-2) 


8.504( 


-6) 


1.243 


-5) 


9.189( 


-6) 





9067 





9154 


-5 


005 


-3) 





7255 


2.495 


1 


453 


1 


154 


-1) 





2695 


6 


962(-2) 


1.532( 


-5) 


2.005 


-5) 


1.304( 


-5) 





8819 





8929 


-6 


999 


-3) 





6432 


2.433 


1 


392 


1 


195 


-1) 





2694 


6 


918(-2) 


1.200( 


-5) 


1.575 


-5) 


7.702( 


-6) 





8616 





8742 


-8 


502 


-3) 





5623 


2.403 


1 


353 


1 


216 


-1) 





2694 


6 


896(-2) 


5.376( 


-6) 


1.547 


-6) 


1.295( 


-5) 





8453 





8591 


-9 


722 


-3) 





4950 


2.372 


1 


305 


1 


238( 


-1) 





2694 


6 


874(-2) 


6.626( 


-5) 


6.375 


-5) 


8.573( 


-5) 





8227 





8382 


-1 


187 


-2) 





4140 




















M/R = 


0.16 vs 


0.16, 


Mb 


= 


1600 vs. 


1600 




















4.263 


2 


858 


5 


643( 


-2) 





2934 


9 


499(-2) 


2.404( 


-5) 


4.283 


-5) 


4.043( 


-5) 





9863 





9903 


-3 


097 


-4) 





9689 


3.653 


2 


452 


7 


007( 


-2) 





2931 


8 


988(-2) 


2.093( 


-5) 


3.718 


-5) 


3.456( 


-5) 





9786 





9833 


-6 


159 


-4) 





9468 


3.044 


2 


037 


9 


028( 


-2) 





2927 


8 


461 (-2) 


1.544( 


-5) 


2.709 


-5) 


2.441 ( 


-5) 





9621 





9679 


-1 


548 


-3) 





8975 


2.739 


1 


821 


1 


044( 


-1) 





2924 


8 


194(-2) 


1.117( 


-5) 


1.921 


-5) 


1.659( 


-5) 





9454 





9522 


-2 


795 


-3) 





8464 


2.556 


1 


684 


1 


147 


-1) 





2922 


8 


036(-2) 


1.693( 


-5) 


1.689 


-5) 


1.253( 


-5) 





9291 





9368 


-4 


126 


-3) 





7937 


2.434 


1 


585 


1 


227 


-1) 





2921 


7 


933(-2) 


1.539( 


-5) 


2.183 


-5) 


1.667( 


-5) 





9127 





9213 


-5 


743 


-3) 





7351 


2.312 


1 


478 


1 


316 


-1) 





2920 


7 


833(-2) 


1.866( 


-5) 


2.815 


-5) 


1.727( 


-5) 





8876 





8986 


-8 


268 


-3) 





6512 


2.251 


1 


405 


1 


365( 


-1) 





2919 


7 


783(-2) 


1.504( 


-5) 


1.274 


-5) 


2.947( 


-5) 





8614 





8745 


-1 


052 


-2) 





5389 


2.220 


1 


355 


1 


390( 


-1) 





2919 


7 


754(-2) 


1.168( 


-4) 


1.423 


-4) 


1.751( 


-4) 





8386 





8535 


-1 


314 


-2) 





4546 




















M/R = 


0.18 vs 


0.18, 


Mb 


= 


1710 vs. 


1710 




















3.654 


2 


692 


7 


173 


-2) 





3106 


9 


902 (-2) 


3.784( 


-5) 


6.963 


-5) 


6.598( 


-5) 





9842 





9887 


-6 


028 


-4) 





9618 


3.044 


2 


247 


9 


228( 


-2) 





3101 


9 


334(-2) 


3.055( 


-5) 


5.616 


-5) 


5.218( 


-5) 





9728 





9781 


-1 


402 


-3) 





9269 


2.740 


2 


018 


1 


066 ( 


-1) 





3099 


9 


046(-2) 


2.478( 


-») 


4.546 


-5) 


4.136( 


-5) 





9617 





9677 


-2 


408 


-3) 





8923 


2.557 


1 


877 


1 


170 


-1) 





3097 


8 


874(-2) 


2.091( 


-5) 


3.724 


-5) 


3.304( 


-5) 





9513 





9579 


-3 


502 


-3) 





8597 


2.313 


1 


678 


1 


340( 


-1) 





3094 


8 


650(-2) 


2.473( 


-5) 


4.064 


-5) 


3.310( 


-5) 





9287 





9367 


-6 


232 


-3) 





7855 


2.191 


1 


570 


1 


441 


-1) 





3092 


8 


542(-2) 


3.136( 


-5) 


5.104 


-5) 


3.558( 


-5) 





9091 





9188 


-8 


975 


-3) 





7201 


2.130 


1 


508 


1 


497 


-1) 





3092 


8 


490(-2) 


2.370( 


-5) 


3.947 


-5) 


2.302( 


-5) 





8934 





9044 


-1 


088 


-2) 





6571 


2.099 


1 


466 


1 


525( 


-1) 





3091 


8 


461(-2) 


1.552( 


-5) 


1.699 


-5) 


3.932( 


-5) 





8783 





8905 


-1 


254 


-2) 





5868 


2.068 


1 


411 


1 


553( 


-1) 





3091 


8 


421(-2) 


1.832( 


-4) 


2.664 


-4) 


3.122( 


-4) 





8547 





8691 


-1 


612 


-2) 





4960 
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TABLE IV: Same as Table [II but for different mass binary systems. The superscript 
massive star. 


" ' " denotes the values for the 


more 


Different mass stars, Irrotational case 
da d n M J \VE{M)\ \VE(GB)\ \VE(FUS)\ 


02/01 03/01 5e c 

II' 111 £ f 

o 2 /Oi a3/ ffl i de c 


X 

x' 















M/R = 


0.12 vs 


0.14, 


M B 


= 0.1299 vs. 


1299 




















6.090 


3.602 


3.199( 


-2) 


0.2572 


8 


790(-2) 


1.249( 


-5) 


2.119 


-5) 


2.011( 


-5) 





9904 





9938 


-1 


076 


-4) 





9819 































9940 





9964 


-5 


955 


-5) 





9889 


3.653 


2.156 


6.637( 


-2) 


0.2565 


7 


191(-2) 


8.519( 


-6) 


1.408 


-5) 


1.267( 


-5) 





9586 





9645 


-1 


116 


-3) 





8950 































9748 





9789 


-5 


560 


-4) 





9355 


3.044 


1.769 


8.579 


-2) 


0.2562 


6 


751 (-2) 


5.310( 


-6) 


8.355 


-6) 


7.043( 


-6) 





9212 





9295 


-3 


284 


-3) 





7858 































9538 





9591 


-1 


541 


-3) 





8732 


2.861 


1.643 


9.359 


-2) 


0.2561 


6 


619(-2) 


8.448( 


-6) 


1.253 


-5) 


9.948( 


-6) 





8985 





9081 


-4 


941 


-3) 





7131 































9420 





9480 


-2 


256 


-3) 





8374 


2.739 


1.553 


9.949 


-2) 


0.2560 


6 


533(-2) 


1.412( 


-5) 


1.903 


-5) 


1.371 ( 


-5) 





8750 





8867 


-6 


731 


-3) 





6375 































9311 





9379 


-2 


988 


-3) 





8038 


2.677 


1.503 


1.027 


-1) 


0.2559 


6 


490(-2) 


1.371( 


-5) 


1.827 


-5) 


1.285( 


-5) 





8583 





8712 


-7 


965 


-3) 





5768 































9243 





9316 


-3 


465 


-3) 





7824 


2.616 


1.446 


1.061 


-1) 


0.2559 


6 


448(-2) 


9.790( 


-6) 


1.443 


-5) 


8.245( 


-6) 





8333 





8481 


-9 


742 


-3) 





4796 































9163 





9242 


-4 


057 


-3) 





7563 


2.585 


1.412 


1.079 


-1) 


0.2559 


6 


428(-2) 


6.266( 


-7) 


5.076 


-6) 


3.032( 


-6) 





8146 





8308 


-1 


118 


-2) 





4100 































9115 





9197 


-4 


414 


-3) 





7404 














M/R = 


0.14 vs 


0.16, 


Mb 


= 0.1461 vs. 


1600 




















4.872 


3.125 


4.591 


-2) 


0.2822 


9 


348(-2) 


2.010( 


-5) 


3.534 


-5) 


3.351( 


-5) 





9865 





9907 


-2 


313 


-4) 





9714 































9913 





9944 


-1 


466 


-4) 





9819 


3.653 


2.346 


6.896 


-2) 


0.2817 


8 


401(-2) 


1.587( 


-5) 


2.772 


-5) 


2.558( 


-5) 





9691 





9746 


-8 


582 


-4) 





9225 































9805 





9846 


-5 


065 


-4) 





9509 


3.044 


1.941 


8.893 


-2) 


0.2813 


7 


901 (-2) 


1.123( 


-5) 


1.927 


-5) 


1.713( 


-5) 





9436 





9507 


-2 


314 


-3) 





8477 































9653 





9703 


-1 


292 


-3) 





9052 


2.739 


1.726 


1.029 


-1) 


0.2810 


7 


650(-2) 


1.417( 


-5) 


2.285 


-5) 


1.876( 


-5) 





9161 





9249 


-4 


405 


-3) 





7630 































9498 





9558 


-2 


366 


-3) 





8577 


2.617 


1.634 


1.096 


-1) 


0.2809 


7 


551(-2) 


1.342( 


-5) 


2.021 


-5) 


1.619( 


-5) 





8974 





9072 


-5 


894 


-3) 





6942 































9405 





9470 


-3 


135 


-3) 





8285 


2.495 


1.534 


1.171 


-1) 


0.2808 


7 


454(-2) 


1.867( 


-5) 


2.751 


-5) 


1.896( 


-5) 





8692 





8817 


-8 


357 


-3) 





6020 































9278 





9353 


-4 


218 


-3) 





7882 


2.434 


1.472 


1.211 


-1) 


0.2808 


7 


406(-2) 


6.809( 


-6) 


1.456 


-5) 


4.008( 


-6) 





8423 





8569 


-1 


046 


-2) 





4925 































9195 





9276 


-4 


967 


-3) 





7611 


2.403 


1.434 


1.232 


-1) 


0.2808 


7 


382(-2) 


2.180( 


-5) 


1.778 


-5) 


3.298( 


-5) 





8210 





8372 


-1 


249 


-2) 





4165 































9146 





9230 


-5 


427 


-3) 





7442 














M/R = 


0.16 vs 


0.18, 


Mb 


= 0.1600 vs. 


1710 




















4.263 


2.988 


5.713 


-2) 


0.3022 


9 


963(-2) 


3.261( 


-5) 


5.931 


-5) 


5.653( 


-5) 





9853 





9897 


-3 


607 


-4) 





9671 































9904 





9939 


-2 


793 


-4) 





9789 


3.044 


2.136 


9.130 


-2) 


0.3014 


8 


885(-2) 


2.260( 


-5) 


4.105 


-5) 


3.774( 


-5) 





9597 





9660 


-1 


765 


-3) 





8917 































9745 





9793 


-1 


225 


-3) 





9308 


2.618 


1.821 


1.123 


-1) 


0.3010 


8 


498(-2) 


1.484( 


-5) 


2.691 


-5) 


2.391( 


-5) 





9312 





9393 


-4 


138 


-3) 





8036 































9578 





9636 


-2 


732 


-3) 





8786 


2.435 


1.676 


1.239 


-1) 
*-} 


0.3008 


8 


336(-2) 


3.135( 


-51 


5.044 


-5) 


3.855( 


-5) 





9072 





9171 


-6 


569 


-3) 





7268 































9448 





9516 


-4 


189 


-3) 





8376 


2.312 


1.568 


1.328 


-1) 


0.3007 


8 


230(-2) 


2.650( 


-5) 


4.263 


-5) 


2.981( 


-5) 





8791 





8912 


-9 


320 


-3) 





6225 































9323 





9400 


-5 


731 


-3) 





7972 


2.282 


1.536 


1.352 


-1) 


0.3007 


8 


204(-2) 


1.853( 


-5) 


3.204 


-5) 


1.787( 


-5) 





8661 





8793 


-1 


047 


-2) 





5663 































9284 





9363 


-6 


185 


-3) 





7842 


2.251 


1.498 


1.377 


-1) 


0.3006 


8 


177(-2) 


1.672( 


-5) 


1.483 


-5) 


3.366( 


-5) 





8461 





8609 


-1 


245 


-2) 





4852 































9240 





9322 


-6 


763 


-3) 





7692 


2.233 


1.473 


1.392 


-1) 


0.3006 


8 


159(-2) 


6.054( 


-5) 


7.237 


-5) 


9.823( 


-5) 





8318 





8478 


-1 


447 


-2) 





4392 































9210 





9294 


-7 


143 


-3) 





7589 



